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The neutral current elastic scattering of neutrinos on Carbon and CH2 targets is 
\Q • computed using the relativistic distorted-wave impulse approximation with relativis- 

tic optical potential. Results for exclusive and inclusive neutrino reactions on 12 C tar- 
get are presented. We show that the nuclear effects on the shape of four- momentum 
transferred squared distribution da jdQ 2 Q E in neutrino neutral-current and charged- 
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current quasi-elastic scattering are similar. We also calculate flux-averaged neutral 



J> , current elastic differential cross section da/dQq E for neutrino scattering from CH2, 

\0 ■ 

' as well as, the neutral-current to charged-current cross section ratio as functions of 

G\ '. 

Qqe- The- value of axial mass Ma is extracted from a fit of da jdQ 2 q E cross section 
measured in MiniBooNE experiment. The extracted value of Ma = 1.28 ±0.05 GeV 



is consistent within errors with the MiniBooNE result. Additionally, for proton ki- 
netic energies above the Cherenkov threshold, the strange quark contribution to the 



neutral current axial vector form factor at Qq E = 0, As, was extracted from a fit of 
MiniBoone data for vp — > vp to vN — > uN cross section ratio. This value is found 
to be As = -0.11 ±0.36 

PACS numbers: 25.30.-c, 25.30.Bf, 25.30.Pt, 13.15.+g 



I. INTRODUCTION 

Neutrino-nucleon neutral-current elastic (NCE) scattering provides an additional infor- 
mation about of structure of the hadronic weak neutral current (NC) and plays an important 
role in searching for the three active neutrino ^active = {v e , v^, v T } conversion to a sterile neu- 
trino v s : a neutrino which has no coupling to neither W ± nor Z° bosons. 

The weak neutral current of the nucleon may be parametrized in terms of two vector and 
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one axial- vector form factors. An additional induced pseudoscalar form factor is presented, 
but its contribution vanishes in the limit of a zero neutrino mass. In particular, the axial- 
vector form factor may be split into a non-strange and strange contributions. The latter 
one is proportional to the fraction of the nucleon spin carried by the strange quarks [1, 2]. 
Thus the axial-vector form factor is crucial for understanding the role that strange quarks 
play in determining the properties of nucleons. 

In order to investigate how the strange quarks contribute to the observed properties of 
the nucleon various reactions have been proposed: deep inelastic scattering of neutrino or 
polarized charged leptons on proton [3, 4], and parity- violating electron scattering [5, 6]. 
The strange vector form factors were measured in parity- violating electron scattering exper- 
iments [7-11]. A combined analysis of these experiments data points to small strangeness 
of the vector form factors [12]. 

Whereas parity-violating electron scattering is sensitive to the electric and magnetic 
strangeness, neutrino-induced reactions are sensitive to the strange quark contribution As 
to the NC axial- vector form factor. A measurement of v{v)- proton NCE at Brookhaven Na- 
tional Laboratory (BNL E734) [13] suggested a non-zero value of As. However, in Ref. [1] it 
has been shown that the BNL data cannot provide a decisive conclusion about the value of 
As when taking into account uncertainties in the vector strange form-factors. Moreover, this 
result suffers strongly from experimental uncertainties due to difficulties in determination of 
the absolute neutrino flux. 

The measurement of the neutral-to-charged-current (CC) quasi-elastic cross section R = 
NCE/CCQE in neutrino-nucleus scattering was proposed in Ref. [14] to extract information 
on the strange spin of the proton because much of the systematic uncertainty is canceled 
by using the ratio. An important effort in this direction was the MiniBooNE experiment, 
that measured the flux- averaged NCE differential cross section da/dQ 2 as a function of 
four-momentum transferred squared Q 2 and ratio R = NCE/CCQE [15]. The MINERvA 
experiment [16] aims at high precision measurements of neutrino scattering cross sections, 
and would be well-suited to examine the Q 2 evaluation of strangeness form factor in NCE 
scattering. 

Recently, the question about an additional sterile neutrino has drawn a considerable 
interest in the literature [18-20]. The short-baseline neutrino oscillation experiment, LSND, 
at the Los Alamos National Laboratory [17] reported evidence of — » v e oscillation, but 



3 



with a squared mass difference Am 2 LSND that is inconsistent within a three neutrino mass 
model with the two other values extracted from solar, atmospheric, and reactor experiments, 
i.e. Am 2 at + Am 2 sol ^ Am 2 LSND . 

One of the most favorable scenarios that accommodates three independent Am 2 values 
is an addition of a sterile neutrino. Because three active neutrinos couple to Z°, the rate of 
neutrino NC events should be unaffected by the three flavor neutrino oscillation. Conversely, 
an existence of a sterile neutrino adds a possibility of a u active — > v s transition that would 
create a deficit in the rate of NC events. 

However, the SNO experiment [21] made a neutral-current rate measurement and showed 
that the total flux of active neutrino from the Sun agree with expectation from the Standard 
Solar model. The Super-Kamiokande experiment excludes — > u s and favors a pure — > u T 
oscillation in its analysis of atmospheric neutrinos where an admixture of the two possibilities 
is allowed [22, 23]. The MINOS collaboration reported [24] the measurements of neutrino 
NC rates and spectra in an accelerator long baseline neutrino experiment. The rates at the 
near and far detectors are consistent with expectations from decay kinematics and geometry, 
providing new support for the interpretation of muon neutrino disappearance as oscillations 
among the three active neutrinos. So, an additional interest in the neutrino-nucleus NCE 
scattering cross section is that this process plays a key role in a search for the parameter 
space available for u active — > u s oscillations. 

It has been shown in Refs. [15, 25, 26] that in order to measure the strange quark con- 
tribution to the nucleon spin using a neutrino-nucleon NCE cross section it is necessary 
to distinguish up — > up from un — > un interactions. Otherwise, the total NCE cross-section 
on both proton and neutron (uN — > z/iV)has a negligible dependence on the nucleon spin's 
strangeness. A detailed analysis of the NCE scattering cross section's sensitivity to the 
strange content of nucleon neutral current was carried out in a relativistic plane-wave im- 
pulse approximation in Ref. [27]. 

Analysis of nuclear structure effects on the determination of the strange quark contribu- 
tion in neutrino-nucleus NCE scattering were performed in Refs. [28, 29] where the relativistic 
Fermi Gas model (RFGM) and relativistic shell model including final state interaction (FSI) 
of outgoing nucleon were used. The effects of FSI were also studied in Ref. [30] within the 
RFGM and in Refs. [1, 31] in the framework of the Random Phase Approximation theory. 
The effects of FSI on the ratio of proton-to-neutron cross section in NCE scattering were 
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discussed in Refs.[28, 32, 33]. 

The effects of FSI on NCE scattering cross section were studied in Refs. [34-36] within 
the framework of a relativistic distorted-wave impulse approximation (RDWIA) with a rela- 
tivistic optical potential. In Refs. [34, 35] important FSI effects arise from the use of optical 
potential within a relativistic Green's function approach. An analys of the sensitivity of NCE 
scattering cross section to the strangeness contribution was presented also in Refs. [37, 38] 
within the RDWIA and relativistic multiple-scattering Glauber approximation. 

In this paper we present the RDWIA calculation of the neutrino-nucleon NCE scatter- 
ing cross section on Carbon and CH 2 . In this approach that was successfully applied in 
Refs. [39-42] to CC quasi-elastic scattering we calculated the flux- averaged da/dQ 2 cross 
section and ratio R(NCE/CCQE) and compare the results with the MiniBooNE data [15]. 
Additionally, the ratio of the predicted event rates in the MiniBooNE high energy up — > nup 
and uN — > nuN event samples was calculated. This ratio is sensitive to the strange quark 
contribution to the nucleon spin, As. Using the MiniBooNE data for this distribution we 
performed a measurement of As and compared it to the MiniBooNE result reported in 
Ref. [15]. 

The outline of this article is as follows: In Sec. II we present briefly the formalism for the 
NCE scattering process and the RDWIA approach. The results are presented and discussed 
in Sec. III. Our conclusions are summarized in Sec. IV. 

II. THE FORMALISM AND MODEL FOR THE NEUTRAL-CURRENT 

ELASTIC SCATTERING 

In this section we consider the formalism for description of NCE exclusive 

u(ki) + A(p A ) -+ u(k f ) + N( Px ) + B(p B ), (1) 

and inclusive 

u(ki) + A(p A ) ^ u(k f ) + X (2) 

scattering off nuclei in the one-Z°-boson exchange approximation. Here k{ = (e,, fcj) and 
kf = (ef,kf) are the initial and final lepton momenta, p A = (sa,Pa)i an d Pb = ( £ b,Pb) 
are the initial and final target momenta, p x = (e x ,p x ) is the ejectile nucleon momentum, 
q = (uu, q) is the momentum transfer carried by the virtual Z°-boson, and Q 2 = —q 2 = q 2 —u 2 
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is the Z°-boson virtuality. As the basic outline follows closely the CC formalism developed 
in Ref. [39], we present a brief review that focuses on those modifications that arise from 
the weak neutral current. 



A. Neutrino-nucleus NCE scattering cross sections 

In the laboratory frame, the differential cross section for the exclusive (anti-)neutrino 
NCE scattering, in which only a single discrete state or narrow resonance of the target is 
excited, can be written as 

R \p x \e x \ Kf \ ^_ L y WfF M ) (3) 



defdVLfdVL x (2ir) 5 £j 
where Qf is the solid angle for the lepton momentum, Q x is the solid angle for the ejectile 
nucleon momentum, R is the recoil factor, G ~ 1.16639 x 10~ n MeV -2 is the Fermi constant, 
Lffi and wjffi are NC lepton and nuclear tensors, respectively. 
The energy e x is the solution to the equation 

e x + e B - m A - u = 0, (4) 

where e B = \/ m B + p\^ Pb = Q ~ Pxi Px = V £ x ~ m2 i an< ^ m A-, m B, and m are masses of 
the target, recoil nucleus and nucleon, respectively. The missing momentum p m and missing 
energy e m are defined by 

P m = Px - Q (5a) 
e m = m + m B — tra (5b) 

From Eq.(4) the total energy of the ejected nucleon is given by 

e x = u + m A -e B ttuj + m-(e m -pl l /2m B ) (6) 

and the nucleon kinetic energy can be written as 

T N = u - (e m - p m /2m B ) w u - e m , (7) 

if one neglects the recoil nucleon energy p 2 m /2m B . As the outgoing neutrino is undetected 
the differential cross section Eq. (3) can be rewritten in "no- recoil" approximation as follows 



dT^dVLfdVLx dsfdQfdQ x 



(8) 
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The leptonic tensor is separated into symmetric and antisymmetric components 

that are given as in Ref. [39]. Note, that the weak lepton NC is conserved for massless 
neutrino and q^Lffi = L^q u = 0. All the nuclear structure information and FSI effects 
are contained in the weak NC nuclear tensor wju, c \ which is given by the bilinear product of 
the transition matrix elements of the nuclear NC operator between the initial nucleus 
state \A) and the final state \Bf) as 

w^ c) = Y.( B f^\ j{ rM)(M nc)] \B f ,p x ), (9) 
/ 

where the sum is taken over undetected states. This tensor is an extremely complicated 
object as, in principle, the exact form for many body wave functions and operators must be 
used. A general mo del- independent covariant form of wjffi and the result of its contraction 
with the leptonic tensor were obtained in Ref. [27] . There it was shown that the contraction 
L^W^^ and therefore the differential cross section in Eq. (3) is completely determined 
by a set of eight structure functions. 

General expressions for the cross sections of the exclusive and inclusive CCQE neutrino 
scattering off nucleus are given in Ref. [39] in terms of weak response functions. In order to 
apply these expressions for calculation of neutrino-nucleus NCE scattering cross sections it 
is necessary to replace G 2 cos 2 9 C — > G 2 , express the response functions as suitable combi- 
nations of the hadron tensor components wfuf\ and calculate the coefficient Vi for massless 
neutrino. The single differential cross section as a function of the outgoing nucleon's kinetic 
energy T N can be obtained after performing integration of the cross section in Eq.(8) over 
solid angles of the outgoing neutrino and nucleon. 

B. Model 

We describe the neutrino-nucleon NCE scattering in the impulse approximation, assuming 
that the incoming neutrino interacts with only one nucleon, which is subsequently emitted, 
while the remaining (A-l) nucleons in the target are spectators. When the nuclear current 
is written as the sum of single-nucleon currents, the nuclear matrix element in Eq.(9) takes 
the form 

(p,B\J^\A) = [ d 3 r exp(it-r)^ (_) (p,r)r^" c ^(r), (10) 
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where r M(nc ) is the NC vertex function, t = e B q/W is the recoil-corrected momentum 
transfer, W = y/ (m A + u) 2 — q 2 is the invariant mass, $ and are the relativistic 

bound-state and outgoing wave functions. 

The single-nucleon charged current has a V—A structure j( nc )^ = Jy nc ^ + J^ nc \ For 
a free-nucleon vertex function, T^ nc "> = Yy nc ^ + T^ nc \ we use the vector current vertex 
function 

r Mnc) = F (nc) {Q 2 )r + ^^F^Q?)/^, (11) 

and the axial current vertex function 

r Mnc) = F K) (g2)r75 + F (-)(Q V75. (12) 

The vector form factors F t (nc) (i = V, M) are related to the corresponding electromagnetic 
ones for proton Ff and neutron F™, plus a possible isoscalar strange-quark contribution Ff, 
i.e. [2] 

4" c) = r 3 (0.5 - sin 2 w )(Ff - F?) - sin 2 w (^f + F?) - F*/2 (13a) 
i> c) = r 3 (0.5 - sin 2 w )(Ff - F 2 ") - sin 2 W (*2 + F«) - (13b) 

where r 3 = +(—1) for proton (neutron) knockout and 9y/ is the Weinberg angle (sin 2 9yy ~ 
0.2313). The axial F A nc ^ form factor is expressed as 

Tf lc) = (r 3 F A - FD/2, (14) 

where F| describes possible strange-quark contributions. In this work we neglect the 
strangeness contributions, i.e., it is supposed that F v = F^ — F A — 0. For the nucleon form 
factors Ff^ the approximation of Ref. [43] is used. Because the bound nucleons are the 
off-shell we employ the de Forest prescription [44] and Coulomb gauge for off-shell vector 
current vertex T v . The vector-axial form factor is parametrized as a dipole with the axial 
nucleon mass M A , which controls the Q 2 dependence of F A . 

The independent particle shell model (IPSM) is assumed in the calculations of the nuclear 
structure, taking into account the short-range nucleon-nucleon (NN) correlation in the 
ground state. According to the experimental data [45, 46] the occupancy of the IPSM 
orbitals of 12 C equals on average 89%. We assume that the missing strength (11%) can be 
attributed to the NN correlations, leading to the appearance of the high-momentum and 
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high-energy component in the nucleon distribution in the target. To estimate this effect in 
the inclusive cross sections, we consider a phenomenological model that incorporates both 
the single-particle nature of the nucleon spectrum at low energy (IPSM orbitals) and the 
high-energy and high-momentum components due to NN correlations. 

For 12 C we use the same relativistic wave functions of the bound nucleon states $ as 
in Refs.[41, 42]. The wave functions were obtained from Ref. [47] as the self-consistent 
solutions of the relativistic Hartree-Bogoliubov equations, derived within a relativistic mean 
field approach. The normalization factors S(a) relative to the full occupancy of the IPSM 
orbitals of 12 C [45, 46] are: S , (lp 3 / 2 )=84%, S(ls 1/2 ) = 100% with an average factor of about 
89%. 

In order to take into account FSI effects in the RDWIA, the distorted-wave function, 
is evaluated as a solution of the Dirac equation containing a phenomenological relativistic 
optical potential. The channel coupling in the FSI [48] of the N + B system is taken into 
account. The relativistic optical potential consists of a real part which describes rescattering 
of the ejected nucleon and an imaginary part that accounts for its absorption into unobserved 
channels. We use the LEA program [49] for the numerical calculation of the distorted 
wave functions with the EDAD1 parametrization [50] of the relativistic optical potential for 
Carbon. This code, initially designed for computing exclusive electron-nucleus scattering, 
was successfully tested against A(e, e'p) data [45, 51], and adopted for neutrino reactions [39]. 

A complex optical potential with a nonzero imaginary part generally produces an ab- 
sorption of the flux. For the exclusive A(l, I'N) channel this reflects the coupling between 
different open reaction channels. However, for the inclusive reaction, the total flux must be 
conserved. In Refs. [52, 53], it was shown that the inclusive CCQE neutrino cross section of 
the exclusive channel A(l, I'N) calculated with only the real part of the optical potential is 
almost identical when calculated via the Green's function approach [52], in which the FSI 
effects on inclusive reaction A(l, I'X) are treated by means of a complex potential, and the 
total flux is conserved. We calculate the inclusive da/dQ 2 sections with the EDAD1 rela- 
tivistic optical potential in which only the real part is included. The inclusive cross sections 
with FSI effects in the presence of short-range NN correlations were calculated using the 
method proposed in Ref. [39]. 
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III. RESULTS AND DISCUSSION 

A. Neutral Current Elastic differential cross section 

The exclusive reaction u+ A — > u+p+B reaction is a good signal sample of (anti) neutrino 
NCE scattering off nuclei. The measurement of the range of the scattered proton, its angle 
with respect to the direction of the incident neutrino 9 P , and its rate of energy loss allow 
to identify the particle as a proton and determine proton kinetic energy T p . In impulse 
approximation, assuming the target nucleon to be at rest inside nucleus, these measured 
quantities (T p , cos 9 P ) determine the neutrino energy through the kinematic relation 

£u = co S e p (i + 2m p /T p y/ 2 -l' (15) 

where m p is the proton mass. In neutrino oscillation experiments with two detectors the 
spectra of the protons as functions of neutrino energies, measured at near and far detectors 
can be used to search for u active — > v s transition that would create a deficit in the rate 
of one proton events at far detector. Note, that precise measurement of the (anti) neutrino 
NCE scattering off neutron appear problematic due to the difficulties associated with neutron 
detection. In Figs.l, 2, and 3 the (anti)neutrino NCE exclusive cross sections op = da p / dT p 
(per bound proton), a n = da n /dT n (per bound neutron), and sum a p + a n are displayed as 
functions of the emitted nucleon kinetic energy for proton, neutron, and proton or neutron 
knockout, respectively. The calculations correspond to Carbon target and incoming energies 
of 500 MeV and 1000 MeV. The upper (lower) panels show the cross sections for neutrino 
(antineutrino) NCE scattering in comparison with results obtained in Refs.[34, 37]. These 
cross sections were also calculated in the RDWIA approach with dipole approximation of the 
nucleon form factors, EDAD1 parametrization of relativistic optical potential and neglecting 
the NN correlations in the ground state of Carbon. We observe that our calculations and 
those performed in the framework of RDWIA formalism are in a good agreement. For 
neutrino (antineutrio) the o p ja n NCE cross section ratio increases almost linearly with 
nucleon energy: from « 0.7 0.7) for T N ss 20 - 50 MeV up to « 0.82 2) for T N = 
700 MeV. 

To study the nuclear effects on the Q 2 distribution, we calculated (with M A = 1.032 GeV 
and F% = 0) the inclusive cross sections (da /dQ 2 ) nuc (per bound nucleon) of the neutrino 
NCE scattering on Carbon. The results for neutrino energies e v = 0.5, 0.7, 1.2 and 2.5 GeV 
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FIG. 1: (Color online) Differential cross sections a p of neutrino (upper panel) and antineutrino 
(lower panel) NCE scattering as a function of the outgoing proton kinetic energy for two the values 
of incoming (anti)neutrino energy: e v = 500 and 1000 MeV, calculated in the RDWIA. The solid 
lines represent the results obtained in this work and the dot-dashed lines are the results of Ref. [34]. 

are shown in Fig. 4 in comparison with cross section for neutrino NCE scattering on a free 
nucleon (da / dQ 2 ) f ree = 0.5[(dcr/dQ 2 ) p + (da/dQ 2 ) n ], where (da/dQ 2 ) p and (da/dQ 2 ) n are 
the cross sections for neutrino NCE scattering on free proton and neutron, respectively. 
Nuclear effects on the shape of the Q 2 distribution, i.e., ratios R(e u ,Q 2 ) = 
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FIG. 2: (Color online) The same as Fig.l, but for neutron knockout (a n ). 

(da/dQ 2 ) nuc /(da/dQ 2 )f ree are presented in Fig.5. The results obtained for neutrino en- 
ergies e v = 0.5, 0.7, 1.2 and 2.5 GeV are compared with those calculated for neutrino CCQE 
scattering in Ref. [41]. We observe that nuclear effects in neutrino NCE and CCQE scatter- 
ing, in general, are similar. The nuclear effects are seen at low Q 2 ; the tail of the momentum 
distribution at high Q 2 , an overall suppression, and slight change in slope in the middle 
region at e u > 1 GeV is also observed. The range of Q 2 where R ~ 1 (i.e., nuclear effects 
are small and therefore cannot affect the measurement of the effective Ma) increases with 
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FIG. 3: (Color online) The same as Fig.l, but for proton or neutron knockout (a p + a n ). The 
dashed lines show the RDWIA results of Ref . [37] . 

the incoming neutrino energy. 

The measurement of the neutral-to-charged-current cross sections ratio in the neutrino- 
nucleus scattering was proposed in Ref. [14] to extract a possible strange-quark contribu- 
tion. Our RDWIA results for R = NC/CC = (da / dQ 2 ) NC / (da / dQ 2 ) cc ratio, obtained 
with Ma = 1.032 GeV and F| = are presented in Fig. 6 as functions of Q 2 for neutrino 
energies e = 0.5,0.6, 1.2 and 2.5 GeV. The inclusive CCQE cross sections (da / dQ 2 ) cc were 
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FIG. 4: (Color online) Inclusive NCE cross sections vs. the four-momentum transfer Q 2 for neutrino 
scattering off 12 C (solid line) and free nucleon (dashed line) and for the four values of incoming 
neutrino energy: e = 0.5,0.7, 1.2 and 2.5 GeV. 

in Ref. [41]. The NC/CC ratio decreases as Q 2 increases from ~ 1.9 at Q 2 ps 0.1 (GeV/c) 2 
and reaches the minimum at large value of Q 2 . The fact that the CCQE cross section 
goes to zero more rapidly than the corresponding NCE one (due to the muon mass) causes 
the enhancement of the ratio at large value of Q 2 close to the upper border of the allowed 
kinematic range of Q 2 . The results obtained in Refs.[34, 37] show similar features. 
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FIG. 5: (Color online) Ratio R(e u ,Q 2 ) vs. the four-momentum transfer Q 2 for neutrino scattering 
off 12 C and for the four values of incoming neutrino energy: e = 0.5, 0.7, 1.2 and 2.5 GeV. As shown 
in the key, the ratios were calculated for neutrino NCE and CCQE scattering. 

B. MiniBooNE flux-averaged differential cross section: comparison with data 

The MiniBooNE collaboration reported [15] high-statistic measurement of the flux- 
averaged NCE differential cross section for neutrino scattering on CH^ as a function of 
Qoe- m this experiment the sum of kinetic energies of all final state nucleons that are pro- 
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FIG. 6: (Color online) Ratio of neutral-to-charged-current cross sections R = NC/CC vs. the 
four- momentum transfer Q 2 for neutrino scattering off 12 C and for the four values of incoming 
neutrino energy: e = 0.5,0.7, 1.2 and 2.5 GeV. 

duced in the interaction T = Jj^Tj was measured and spectrum of NCE events dN^cE/dT 
was reconstructed as a function T ~ u. Assuming that the target nucleon is at rest, Qq E 
was determined for each event as 



Q 2 QE = 2mT = 2m^2 T { . (16) 
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FIG. 7: (Color online) Flux-averaged da jdQ 2 Q E cross section per nucleoli (upper panel) for neutrino 
scattering on CH2 and NCE/CCQE cross section ratio (lower panel) as a function of Qqe- The 
NCE cross section and CCQE/NCE ratio calculated with values of Ma = 1-28 GeV (blue line), 
1.37 GeV (red line), both with the value As = 0.0. Also shown is the strange quark effect on the 
NCE cross-section and the ratio with a value of As = —0.11 and Ma = 1-28 GeV (black line). The 
MiniBooNE data are shown as points. 
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MiniBooNE reported the NCE differential cross section in the range of Qq E from 0.1 to 
about 1.65 (GeV/c) 2 . This differential cross section distribution was fitted with a dipole 
axial form factor and best fit for M A = 1.39 ± 0.11 GeV was obtained. Using the data from 
charged- current neutrino interaction sample, the NCE/CCQE cross sections ratio as a 
function of Qq E was measured. One should understand that, in fact, the NCE flux- averaged 
differential cross section da/dT was measured, which under the assumption of a scattering 
on free nucleon, was recalculated as da/dQ 2 QE = (da/dT) /2m. The details of the unfolding 
procedure in measurement of the cross section can be found in Ref. [25]. 

We calculated the NCE flux-averaged inclusive (da NC E/dT) cross section in the frame- 
work of RDWIA approach that was recalculated as (da^cEl 'dQ 2 Q E ) = (dfNCE/dT) /2m. The 
details of the calculation (da NC El 'dT) cross section are described in Appendix A. To extract 
a value for the parameter Ma we calculated this cross section with the Booster Neutrino 
Beamline flux [54] using the Q 2 q E bins AQ 2 = Qf +1 — Q 2 similar to [15] 

and NCE to CCQE cross section ratio 

R t = NCE/CCQE = {da NCE /dQ 2 QE ) t / (da cc /dQ 2 ) v (18) 

The CCQE differential cross section dacc/dQ 2 was calculated in the RDWIA approach 
in Ref. [42]. The fit to the extracted flux-averaged {da NC El ^Qqe) yield the parameter 
M A = 1.28 ± 0.5 GeV. Fig. 7 shows the MiniBooNE measured flux-averaged differential 
daNCEl t dQ 2 q E cross section and R = NCE/CCQE ratio [15] as a function of Qq E compared 
with the RDWIA calculations with the value of Ma = 1.28 GeV. The result obtained with 
Ma = 1.37 GeV, that was extracted in Ref. [42] from the fit to measured in Ref. [55] 
flux-integrated CCQE cross section dacc/dQ 2 , also is shown. 

There is an overall agreement within errors between the RDWIA predictions and the 
MiniBooNE data at high Q 2 QE . However one should note that at Q 2 QE > 0.4 (GeV/c) 2 the 
inclusive cross section, as well as, R = NCE/CCQE ratio, calculated with the value of 
Ma = 1-28 GeV has a better agreement with data. At Q 2 q E < 0.25 (GeV/c) 2 the calculation 
underestimate both cross section and the NCE/CCQE ratio by 17% or less. 
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C. MiniBooNE up —y up/uN — ^ vN differential cross section cross-section ratio: 

comparison with data 

In addition to the uN — y uN differential cross-section, MiniBooNE has published the 
up — y up to uN —y uN ratio at Q 2 > 0.7 GeV 2 (above the Cherenkov threshold for protons 
in mineral oil) [15]. This result is interesting, because it should be sensitive to the strange 
quark contribution to the axial form factor. MiniBooNE has reported the measurement 
As = 0.08 ± 0.30, based on a Nuance prediction [56]. 

The up —y up/uN —y uN ratio was reported as a function of the MiniBooNE recon- 
structed nucleon kinetic energy T rec . Also the migration matrices were published in Ref. [57], 
which carry the detector resolution and efficiency information. Using them one can smear 
the predicted cross-sections and obtain the predicted event rates as a function of T rec . The 
procedure for carying out calculations of event rates in terms of the MiniBooNE recon- 
structed energy is described in an Appendix of Ref. [25]. We performed the calculation of 
the up —y up/uN —y uN ratio based on our neutrino interaction model and compared it 
to the MiniBooNE data. We have calculated our prediction of the event rates for different 
values of As covering the range from —0.4 to 0.4. An example of the calculation is shown 
in Fig.8. 

Using the full error matrix for the ratio published in Ref. [57] we calculated the x 2 distri- 
bution between data and the MC. Our calculation leads to: 

As = -0.11 ±0.36. 

with Xmin = 33.4 for 29 degrees of freedom. This result is consistent with all other mea- 
surements of As, including the one reported by MiniBooNE. We show a calculated NCE 
cross-section and NCE/CCQE ratio with the values of As = -0.11 and M A = 1.28 GeV in 
Fig. 7. As one can see the effect of strange quarks is small, but the agreement between data 
and our prediction does improve a little bit at low Q 2 region with As = —0.11. 

IV. CONCLUSIONS 

In this article we study neutral-current elastic (anti) neutrino scattering on Carbon and 
CH 2 targets in the framework of RDWIA approach placing particular emphasis on nuclear 
effects. We calculated the NCE exclusive da/dT cross sections for nucleon knockout in 
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FIG. 8: (Color online). MiniBooNE up — > up/uN — > vN ratio as a function of T rec . The prediction 
for As = 0.0 and As = 0.4 are shown as red solid and blue dashed histograms, respectively. The 
MiniBooNE data are shown as points with the total error bars on top of them. 

(anti)neutrino scattering on a Carbon target. The calculation presented in this paper are 
consistent with the RDWIA cross sections of Ref. [34, 37]. 

We also calculated the da/dQ 2 inclusive cross sections for neutrino scattering on 12 C, as 
well as, on free nucleon for different neutrino energies and estimated the range of Q 2 where 
nuclear effects on the shape of Q 2 distribution are negligible. We show that these effects in 
the CCQE and NCE scattering are similar. 

Using the RDWIA approach with the Booster Neutrino Beamline flux [54] we extracted 
axial mass from a "shape-only" fit of the measured flux-averaged do jdQ 2 Q E differential cross 
section. The extracted value of Ma = 1-28 ± 0.05 is in agreement within errors with the 
MiniBooNE result of Ma = 1.39±0.11 GeV. There is a good overall agreement within errors 
in the range of 0.25 (GeV/c) 2 < Q 2 < 1.65 (GeV/c) 2 between RDWIA prediction and the 
MiniBooNE data: the measured MiniBooNE NCE flux-averaged differential cross section 
da/dQq E on CH 2 and the NCE /CCQE cross section ratio. However, in the range of low 
Q 2 < 0.25 (GeV/c) 2 the calculations underestimate the measurements by 17% or less. 

Using MiniBooNE data for the high energy up — > up to uN — > uN ratio the value of 
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As = —0.11 ± 0.36 has been extracted based on our model, which is consistent with other 
measurements of As. 

We conclude that the RDWIA approach was successfully tested against neutrino CCQE 
and NCE scattering on 12 C. 
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Appendix A: Flux-averaged NCE inclusive cross section 

MiniBooNE measured the flux-averaged NCE differential cross section (per nucleon) on 
CH 2 , averaged over three process: scattering off free protons in Hydrogen, bound protons 
in Carbon, and bound neutrons in Carbon. This cross section can be expressed as 
dcr V N _ \q ^ da PiH ^ 3^ ^ da P:C 3^ ^ da njC 
^Qqe ^ ^Qqe ^ ^Qqe 7 dQq E 
where da P: n is the NCE cross section on free protons (per proton), da P: c is the cross section 

on bound protons (per proton), da n ^ is the cross section on bound neutrons (per neutron), 
Cu P ,h, Cu P ,c, and C m> c are the efficiency correction functions, that are given in Table IV of 
Ref. [15]. 

In this paper the NCE inclusive cross sections for neutrino [y^ + u e ) scattering on bound 
proton da v j dQ 2 Q E and on bound neutron da n j dQ 2 Q E are calculated in the RDWIA approach. 
The flux-averaged (da / dQ 2 Q E ) cross section can be written as 
da 2 \ _ 1 p— da 

?QE J ~ *i,„,„ *%, 

where I Vii (I Ue ) is the neutrino spectrum and $ is the neutrino flux [y^ + u e ) in v- mode of 

beam, integrated over < e v < 2.6 GeV. This definition of the flux-averaged NCE inclusive 

cross section is similar to the definition in Ref. [55] of the flux-integrated CCQE differential 

cross section da/dQ 2 . 



\ 1 /" £max da- 

(I() 2 (QIe) ) = $ / -t^{Qqe^)\K^) + lM]de v , (A2) 

^ QE I J £"min QE 
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